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SUMMARY 

The digital derived rate (DDR) algorithm Is used to calculate the rate of 
rotation of the Centaur upper-stage rocket. The DDR Is a highly nonlinear 
algorithm, and classical linear stability analysis of the spacecraft cannot be 
performed without linearization. The performance of this rate algorithm Is 
characterized by a gain and phase curve that drop off at the same frequency. 
This characteristic Is desirable for many applications. 

This paper deals with the Investigation of a linearization technique for 
the DDR algorithm. The linearization method Is described. Examples of the 
results of the linearization technique are Illustrated, and the effects of 
linearization are described. A linear digital filter may be used as a sub- 
stitute for performing classical linear stability analyses, while the DDR 
Itself may be used In time response analysis. 


INTRODUCTION 

The powered-phase autopilot (ref. 2) Is a digital control system for the 
Centaur vehicle (fig. 1). It Is used to control the vehicle during full thrust 
operation. Control Is maintained by using hydraulic actuators to rotate, or 
glmbal, the main engines. This control system Is Illustrated by block diagram 
for one plane of motion (the pitch plane) In figure 2. 

The attitude error of the Centaur vehicle (the difference between actual 
and desired flight direction) Is sampled and digitized according to the A/D 
conversion In figure 3. The attitude error Is then available In the discrete 
domain. 

The control law Is a proportlonal-lntegral-derl vatlve (PID) controller 
with a first-order digital Integrator. The autopilot operates at a sampling 
frequency of 50 Hz. Discrete rate Is a function of attitude error (n), sam- 
pling period (T), and A/D quantization resolution (RES). 

. RES(,, i - Vi> 

n = j (1 ) 

RFS Is equal to 0.02919 deg/b, so the minimum change In attitude error 
over one sampling period results In a rotation rate of 1.4595 deg/sec. Rate 
changes of <0.1 deg/sec are desired during one period. For this reason, a 
50 Hz autopilot must use a higher-order algorithm to maintain good rate Infor- 
mation during transient conditions (engine startup, etc.). There are also 
limits on available memory storage and execution time in the digital flight 
computer. 


Centaur flights use the DDR algorithm to calculate rotation rate by 
tracking changes In attitude error. The proper choice of DDR constants 
(ADOR, FDOR, TAU) will give good Information at normal vehicle frequencies. 

The block diagram for the DDR algorithm Is shown In figure 4, and Its con- 
stants are defined In the nomenclature. 

This algorithm has been used successfully In previous Centaur missions. 
Analytical performance evaluation on these past missions has been dependent on 
time response simulation. Classical frequency domain analysis requires linear- 
ization. Linear analysis computer codes such as EASY5 (ref. 5) use simple 
single-input, single-output blocks to model elements such as the DDR. EASY5 
could not be used to model the system In figure 2, because no connection 
between Input and output was recognized within the DDR algorithm. 

To perform linear analysis on the powered-phase autopilot, It was neces- 
sary to find a linear approximation to the DOR. The problem was resolved by 
finding the Fourier transform of the Input and output of the DDR algorithm 
for a simple quantized sine-wave Input. The frequency spectra of Input and 
output were used to find a frequency response spectrum. The nonlinearity of 
the DDR algorithm made It easier to use a series of sine waves of different 
frequencies as Input to the algorithm than to use a series of random Inputs to 
find the frequency response. 

The frequency response of the DDR algorithm Is shown In figure 5. It 
Is a function of the Input amplitude, the DDR constants ADDR, FDDR, and 
TAU, and of noise. 

This work will Illustrate: (1) the usefulness of the Fourier transform 

In solving this sort of nonlinear problem, and (2) the criteria used In 
selecting a linear approximation to this type of element for use In analytical 
simulation. 


NOMENCLATURE 

A amplitude of Input sine wave, deg 

ADDR DDR Increment limit, deg/sec 

DDR digital derived rate, deg/sec 

F resolution of FFT output, Hz 

FDDR DDR Increment for direction change, deg/sec 

K DDR Incremental Index 

N number of samples In record 

RDDR reference derived rate, deg/sec 

RES A/D quantization resolution, deg/b 

T sampling period, sec 
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t time duration of record, sec 

TAU DDR multiplicative scale factor 

n current attitude error, b 

n actual rotation rate, deg/sec 

no attitude error, deg 

n|c reference attitude error, b 

u frequency of Input sine wave, rad/sec 

cutoff frequency of linear DDR approximation, rad/sec 
5 a damping ratio of linear DDR approximation 


THE DIGITAL DERIVED RATE ALGORITHM 

The operation of the DDR algorithm Is based on constants FDDR, ADDR, 
and TAU. Attitude error In bits (n) Is Input to the algorithm and digital 
derived rate (DDR or n) In deg/sec Is output. 

FDDR Is a rate Increment used In the case of a change In vehicle direc- 
tion. If the difference In attitude error between current and reference values 
(n - n|() Is negative and the current rate (DDR) Is zero or positive, then 
DDR Is decreased by FDDR. If the error difference Is positive and DDR Is 
zero or negative, then DDR Is Increased by FT'DR . 

ADDR Is a rate Increment limit used when the difference In attitude error 
and DDR are both nonzero and of the same sign. A reference rate Is found by 
equation (2). 


RES(n - nj 

RDDR = — (2) 

The current value of DOR Is not permitted to be Incremented by more than 
ADDR, nor Is It permitted to be a value less than RDDR. This provides a good 
deal of attenuation of noise and higher frequency components, since DDR Is 
not permitted to become too large. 

TAU Is a scale factor used tc keep DDR from becoming too large In the 
case of a change In vehicle direction. If the vehicle did not change direc- 
tion, then the counter (K) Is Increased by one and DDR Is not changed. 
Otherwise, DDR Is multiplied by TAU (which must be <1 ) . 

The counter K defines the number of sampling periods that a value of 
reference attitude error (r^) Is delayed. When K Is Incremented, the 
change In attitude error over K sampling periods has not been significant. 
When the change In attitude error becomes greater than 1 dB, then DDR Is 
changed according to the rate Incremerts, the current attitude error becomes 
the new reference and K Is reset to onr?. 
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FREQUENCY RESPONSE OF THE DOR ALGORITHM 

The steps In solving for the frequency response of the DDR algorithm 
are shown by flowchart In figure 6. 

(1) Input a sine wave of given amplitude, frequency, and time duration to 
the DDR algorithm with quantization. For a given Input sine wave, 


rip = A sin «t 

the Ideal derivative Is given by 

n > uA cos ut 

- wA sin (ut ♦ 90°) 

The transfer function for the Ideal differentiator Is 

n w sin (ut ♦ 90°) 

9(“>) = r - — 


(3) 


(4) 


(5) 


The magnitude and phase of this transfer function are 

|g(u) | = 20 log u(dB) 
Zg(u) * 90° 


( 6 ) 


Since the attitude error Is harmonic, the frequency response of the DDR 
should approximate equations (6) for low frequencies. 

(2) Find the DDR output for the given Input. A sample of DDR time 
response Is Illustrated In figure 7. For this cas?, the DDR constants were 
chosen as follows: 


ADDR = 0.5 deg/sec 

FDDR = 0.2 deg/sec (7) 


TAU = 0.5 


along with resolution 0.02919 deg/b and sampling time 0.02 sec. The noise 
level used was +1 dB (RES). The Input was defined by the following (used In 
eq. (3)) 


A - 0.5° 

u = 1.26 rad/sec 

From this time response plot, the magnitude and phase of the DDR 
function can be read directly with some accuracy. 

|g(u) | = 1.8 dB 

zg(u) = 80° 
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( 8 ) 


transfer 


(9) 
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which agrees fairly well with equations (6). However, choosing frequency 
Information from time responses Is less reliable as noise becomes more dominant 
In the time response. 

(3) Take the FFT of both Input and output signals. The FFT of the 
Input and output signals (1 and 4) will give the frequency spectra of those 
signals. If the sampling frequency Is given, the time duration of the signal 
(t) must be chosen to give a desired resolution (F) of the frequency spectra. 
For this example, F was chosen to be 0.01 Hz, so t was 100 sec. The number 
of points In the record (N) was then taken to be equal to the time duration of 
the signal divided by the period, or, 5000 points. The FFT returns 2500 
valid frequency domain points, and the range of the spectrum was normalized to 
DC to 25 Hz. 

(4) Choose the gain and phase at the Input frequency. From the points 
found In step C, the gain and phase of the DDR are chosen at the point where 
the frequency Is the same as that of the Input sine wave. 

(5) Repeat steps 1 through 4, Increasing the Input frequency of the sine 
wave. This step gives a sweep of frequencies of Interest. Frequencies higher 
than about 7.5 Hz were filtered prior tc sampling, and the response for fre- 
quencies >5 Hz v.as not calculated. Tie frequency response shown In figure 8 
Is for an Input amplitude of 0.5°, *or cases both with and without noise. 

(5) Repeat steps 1 through 5 for a selection of Input amplitudes of 
Interest. For most Centaur flights, attitude error limit cycles are <2°. 
Frequency responses were taken at Input amplitudes of 0.1, 0.3, 0.5, 1, and 2°, 
to get a good representation of this range. 

A FORTRAN program was written to perform the steps listed. The Input sine 
wave was sampled and quantized according to the plot In figure 3. In some 
cases, random noise was added to the quantized signal prior to calculation of 
DOR. Some cases were also run without noise. Time responses were generated 
with another program to check the results. 


RESULTS 

The frequency responses shown In figure 8 are typical for many combina- 
tions of DDR parameters and Input amplitudes. For all cases observed, the 
gain and phase began to drop off at the same frequency. The character of the 
response was also seen to be affected by noise. 

The gain and phase plots both give the Impression that the DDR may be 
modeled as an Ideal differentiator with a second-order lag. 


G(s) = p . (10) 

S . 2f A ^ A S * * A 

If the approximation Is chosen to fit the phase, however, the cutoff fre- 
quency will be about a decade greater than If the approximation Is chosen to 
fit the gain. The approximation Is therefore accurate over a small frequency 
range of Interest. 
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The nonlinearity caused Inconclusive results for fitting an equation (10) 
approximation to the data by Levy's method (ref. 2). 

The DOR algorithm as modeled according to the approximation of 
equation (10), with u/^ set equal to the cutoff frequency observed In 
figure 8, and 5a chosen to be appropriate for the phase characteristic 
for the case with noise. 


«A * 9.42 rad/sec 
*A = 0.7 


( 11 ) 


This approximation Is accu ate to within 1 dB In gain and within 250 to 30° In 
phase In the vicinity of the cutoff frequency. Below the cutoff frequency, the 
accuracy Is even better. For the case of most Centaur flights, where attitude 
error limit cycles oscillate at <1.5 Hz, this Is a good approximation for 
linear analysis. 


CONCLUSIONS 

The method shown here generated a suitable frequency response for a 
linear approximation to the DDR algorithm. The approximation was used In 
simulations written In EASY5 for root locus and frequency domain analysis. 

The method could be used to find frequency response Information for other 
nonlinear components, provided a certain amount of Information Is known about 
the fundamental behavior of the component, so the right parameters may be 
varied, and so the results obtained by this method may have meaning. For 
Instance, It was necessary In this case to know that the effects of the three 
DDR constants were coupled, anJ that the response of the DDR might vary 
with Input amplitude. 

This technique Is only necessary because of the nonlinear nature of the 
DDR algorithm. If a discrete linear algorithm could be developed that would 
maintain accuracy at low frequencies and high noise levels, without exceeding 
computer storage limitations, It would be a great aid to linear stability 
analysis . 

The approximation constants «a and 5 a were found by this technique 
as functions of the DDR constants. Input amplitude, and noise. This would be 
a more powerful analysis tool If DDR constants could be calculated as 
functions of «a and 5 a, either analytically or graphically. 
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